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SOUND DIFFRACTION AT WALL IMPEDANCE DISCONTINUITIES IN A CIRCULAR 


CYLINDER - INVESTIGATED USIN6 WIENER-HOPF TECHNIQUE 


Y.-C. Cho* 

National Aeronautics and Space Administration 
Lewis Research Center 
Cleveland, Ohio 44135 


ABSTRACT 

Rigorous solutions are presented for sound diffraction In a circular cyl- 
inder with axial discontinuities of the wall admittance (or Impedance). Analy- 
tical expressions are derived for the reflection and the transmission coeffi- 
cients for duct modes. The results are discussed quantitatively In the limits 
of small admittance shifts (a) and of low frequencies (lea). One of the more 
remarkable results Is the low frequency behavior of the reflection coefficient 

( R oo) of the fundamental mode. For the mode of a hardwall duct reflected from 

the junction with a softwall duct, R^ ♦ - (l - vlca V27 T" a); this result Is 

In contrast to the frequency dependence of the reflection from the open end of 

a hardwall duct, for which Rqq ♦ - [l - (ka) 2 /2]. 

NOMENCLATURE 

duct radius 
sound speed 
see eq. (47) 

Green's function; x and Xq, respectively, being coordinates of 
observation and source points 

Bessel function of order m 

u/c, free space wave constant 

propagation constant of (m,&) mode 

see eq. (49) 

see eqs. (55) and (56) 

normalization factor, see eq. (7) 
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p acoustic prossurt 

normalized eigenfunction, see eq. (6) 

rJJj reflection coefficient for (m,n) mode Incident and (m,J) mode 

reflected 

(r,f,x) cylindrical coordinates of observation point 

( r O* , D* x o) cylindrical coordinates of source point 

S surface, see fig. 2 

Tj v (l ♦ 2) transmission coefficient for (m,v) mode of duct 2 excited by 
Incidence of (m,n) mode of duct 1 
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mi 
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transmission coefficient of (m,n) mode of duct 1 excited by 
Incidence of (m,v) mode of duct 2 

time 

particle velocity 

Neumann function (Y - Bessel function) of order m 

eigenvalue corresponding to (m,i) mode 
complex variable, real being conjugate to spatial coordinate 
variable x 

Fourier transform of see eq. (46) 
reduced Green's function, see eq. (28) 

- nj, admittance shift 
Dirac delta function 
admittance (pcu/p) 
see eq. (C6) 
see eq. (8) 
density of medium 
Fourier transforms of \ 
velocity potential 
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$m eqs. (33) tnd (34) 
angular frequency 
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INTRODUCTION 

Lined ducts have long been employed as silencers of noise generated In a 
fan duct system. The sound attenuation In a long, uniformly lined cylinder Is 
adequately understood In terms of acoustic nave modes with complex propagation 
constants. However, In practice the liners are finite In length, and are some- 
times made with segments thus providing discontinuities of the wall Impedance. 
At the Impedance discontinuities, a wave undergoes diffraction, being partly 
reflected and partly transmitted. The reflection accompanies the mode scatter- 
ing, and the transmission Involves the conversion of the mode to a new set of 
modes. In addition to the sound absorption by the liners It has been suspected 
that the reflection and the mode conversion can also contribute appreciably to 
the noise reduction. Such a contribution still requires a comprehensive 
study. .This paper presents a rigorous treatment using the Wlener-Hopf tech- 
nique^ of the diffraction of a mode due to the discontinuity of the wall Im- 
pedance In a circular duct. The analysis Involves an Incident wave of an arbi- 
trary mode, and can be expanded to a case of multiple diffractions which occur 
In a duct Involving more than one axial discontinuity In the wall Impedance. 

Similar approaches were used previously in studies of the diffraction of 
the fundamental mode In rectangular ducts. I 3-5 ) The formal solution Involves 
Infinite product terms, and can be easily manipulated to produce numerical re- 
sults In some limiting cases with the fundamental mode Incident. However, with 
a higher order mode Incident, the numerical evaluation requires calculation of 
complex eigenvalues of many modes, and the Infinite product terms converge 
weakly especially when the Impedance discontinuity Is large. In the present 
analysis the formal solution Is obtained first, and Is used only to determine 
the analytical properties of various functions Involved In the Intermediate 
stages of the computation. 


FORNULATION AND SOLUTION OF WIENER-HOPF TYPE INTEGRAL EQUATION 

Nodal Solution 

This section begins with a review of modal solutions In a circular duct. 

A duct mode corresponds to an elgensolutlon of the steady wave equation In a 
uniform duct which 1$ subject to a homogeneous boundary condition on Its wall. 
With simple harmonic time dependence e~ lMt of the field, the wave equation 
Is reduced to the scalar Helmholtz equation, which Is written for cylindrical 
coordinates (r,*,x) as follows: 

(1 > 

Here k • «/ c, and v Is the acoustic velocity potential and accordingly, the 
particle velocity u and the acoustic pressure p are given by 
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U « V¥, 


p = ipckt, 


( 2 ) 

(3) 


where p is the steady density of the medium. The wall boundary condition is 
specified in terms of the specific acoustic admittance n : 


11 

3r 


r=a 


- r=a ' °- 


(4) 


An eigensolution to this problem is obtained as: 

• ikx 

» - e ,m \,(r)e m . 


(5) 


Here Q m «(r) is the normalized Bessel function: 


~ IT" ^m (°mi a) 
nu v ' 


( 6 ) 


with the normalization factor given as 


ITU 


aJ ( a ) 
m v mr 

v M mi 


(7) 


where 


M m0 = 1 for a = 0, 
mi mi 


1 - 


m 


, 2 + 


(kan)^ 


n-1/2 




for o 4 0. 
mi 


( 8 ) 


The solution in equation (5) corresponds to the (m,i) mode, m and «, being 
integers called the circumferential and radial mode numbers, respectively. The 
constant a is the corresponding eigenvalue that is obtained as the ith 

root of the equation 


60 la) 

a — ^ ikanJ m (a) = 0, (9) 

and k mt is the propagation constant given by 
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For an acoustically absorbing wall for which /?«(n) > 0, the phases of a 

and k are chosen such that 
m 

A*(a m )>0, U{ a mt ) < 0, (11a, b) 


*‘( k mi> 


> 0, 



) > 0. 


(12a, b) 


It follows then from inequality (12a, b) that the solution in equation (5) cor- 
responds to the (m,i) mode propagating to the positive x-direction. For the 
propagation to the negative x-direction, -k m , is used in place of k mi . Note 
that the lowest radial mode number is here chosen to be zero, corresponding to 
the zeroth (or smallest) root of equation (9). 

The problem to be treated here is concerned with the diffraction of a duct 
mode due to the axial discontinuity of the wall boundary condition in a circu- 
lar cylinder, as illustrated in figure 1. The cylinder is, for convenience, 
thought to be composed of two duct elements coupled at x =- 0. The elements 
are designated by 1 and 2, respectively, for x < 0 and x > 0. The respec- 
tive wall admittances are and r^: that is, for x < 0 


3 V 

3r 


r=a 


ikn-it = 0, 
1 r=a 


(13) 


and for x > 0 


3V 

3r 


r=a 


ik^f 


r=a 


= 0. 


(14) 


The eigensolutions are 




> im? (rip 
mu v 


ikil’x 


for x < 0, 


(15) 




ik (2) x 

e im<P Qm! } (r)^ mt 


for x > 0. 


(16) 


Here Q „ and k are superscripted with 1 and 2, respectively for the 

duct elements 1 and 2. The corresponding eigenvalues are denoted by a' 1 ' 

( 2 ) mt 
v ' However, the circumferential mode number m is a conserved quan- 


and 


mi 


tity because the cylindrical symmetry is maintained throughout the duct, and 
the angle dependence e im<p will be often omitted. 


Consider an incident wave which is composed of a single mode of the duct 
element 1, propagating towards the center (x = 0). Upon arrival at x = 0, 
the wave is partly reflected and partly transmitted. With the (m,n) mode in- 
cident, the resultant wave can be written as, for x < 0 
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* = 


ik (D x 

q ( 1 ) (r )e mn + 
y mn ' ' 


-1k (1 >x 

R nj<l>Qmj )(r ) e ^ ’ 


and for x > 0 


oo 

X 

v=0 


t"J (i 

nv 


ik (2) x 
2)Q^ ) ( r >e m . 


( 17 ) 


(18) 


The first term in equation (17) is the incident wave, and the second term is 
the reflected wave, which is composed of many modes of the duct element . 
Equation (18) represents the transmitted wave which is composed of many modes 

of the duct element 2. The constants R™.(1) and T™ (1 ♦ 2) are the reflec- 

nj' nv' 

tion and the transmission coefficients. The first subscript of the coeffi- 
cients is the radial mode number of the incident wave, and the second subscript 
the radial mode number of a mode in the reflected or the transmitted wave. 

These coefficients will be here determined by using the Wirner-Hopf technique. 


Wiener-Hopf Integral Equation 


For the formulation of the Wiener-Hopf type integration, we consider the 
Green's function, G(x,Xg), which is a solution of the equation 


1 i_ r !_ + i_ iL. + iL. + k 2 

^ r3r 3r 7 7? 77 


G(x,x Q ) = -p «(r - r 0 )6(<p - 


6(x - x n ). (19) 


Let it be subjected to the boundary condition of the duct element 1; the 
Green's function is then constructed in terms of the eigensolutions for the 
duct element 1 : (6) 


G(x 


•*o ) * S 


im( <p-<pq) 


m= — o 



r 1k mi ) | x - x ol 
— jn e 

m 


( 20 ) 


In the absence of sources, the acoustic field inside the duct can be 
expressed in the integral form 

¥(x) =f J^(xq) v Q G(x ,x Q ) - G(x,x 0 )v 0 v(x 0 )J • n dS Q . (21) 


The boundary surface S is shown in figure 2 and n is the outward unit vec- 
tor normal to the surface. It is convenient to divide the surface integral 



0F poor 


OOriilfY 


Into three parts: the integrations on the duct wall, Si, on the cross section 

at xq - -«*, S^, and on the cross section at xq « •, S3, as follows: 

» - Ij + I 2 + I 3 . (22) 



The integral 1 3 is zero because t and G decay exponentially as xq ♦ •, 
as one can see from equations (18) and (20). As for the integral I 2 , t In 
equation (17) is used. The contribution from the reflected wave is zero based 
on an argument similar to the argument for I 3 . With the expression for the 
incident wave and equation ( 20 ) inserted into equation (24) one can readily 
obtain 

* e ’ mf Q il (r > e mn • (26) 


This is none other than the incident wave. Now turn to the integral Ij. For 
xq < 0, t and G satisfy the same boundary condition given in equation (13). 
Thus, the integrand is zero for xq < 0, and we have 



where xm ’ s the educed Green's function. 


( 27 ) 
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x-x. 


(28) 


It is convenient to decompose y as 

y = v i + Yq. (29) 

Here is the incident wave. The function Yp represents the reflected 

wave for x < 0, and also the transmitted wave minus the incident wave for 
x > 0 . Note that for x > 0 neither y^ nor Yp satisfies the boundary con- 
dition in equation (14), although y does. The condition for Yp can be 
readily derived as, for x > 0 


= ikr^Y - iknjY^ (30) 


Using equations (29) and (30), and the fact that the incident wave and the 
Green's function satisfy the same boundary condition, one obtains 




>1 ■ - 1ka ^ 

,a;x 

- Xp) Y(a,x 0 

)dx 0 

(31) 

where a = n ^ 

- n,. It follows from equat 

ions 

(22), (26), 

and (31) and 

h * 0 

that 







Yp ( r, 

x) = -ika a / 

•'o 

d V. |r '*i< ■ 

- x o ) 

[r ( (a,x 0 ) * 

»D (a>x 0>]- 

(32) 

Introduce the 

functions y±(x) 

: 






1 

r Y D (a,x) 

for 

x > 0, 




Y + ( X ) = <( 





(33) 



lo 

for 

x < 0. 




| 

r o 

for 

x > 0, 




Y_(x ) =< 





(34) 



*D (a ' x) 

for 

x < 0. 




Then Equation (32) can be written for r = a as 


t + (x) + y_(x) = -ika a f dx 0 Y m (a,a;x - Xq) ( a,x Q ) + y + (x q )]. 


( 35 ) 


This is a Wiener-Hopf type integral equation. 
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Fourier Transform Solution 
First consider the Fourier transforms of t*(x): 

•+(b) - — — f T + (x)e 10x dx, (36) 

vm -- 


♦ (6) - -±- r t (x)e iBX dx. 

vnr -- 


(37) 


With the lowest radial mode number corresponding to the least attenuating mode, 
the field behaves asymptotically for x ♦ as 


mo 


(38) 


Note that, if is less than is replaced by k^ (the 

L mn J L mn J ’ mo mn 

incident wave constant) here and in equations (40) and (42) below. And for 

x ♦ — 


-ik»>x 

mo 


(39) 


It follows that *+(b) am«J *_(b) are analytic, respectively, in the upper half 
plane (UHP) and in the lower half plane (LHP) of the b space. Here the half 
planes are defined as follows: 


U(b) > - A»[k^^] for UHP, (40) 

U(b) < 4*[k^] for LHP. (41) 

Note that the two half planes overlap in the region specified by 

-^[ k io } ] < jW(B) < A *[ k TO ) ]* 

The functions *+(b) and *_(b) are both analytic in this region, and the in- 
verse transforms are obtained as 


t + (x) 
» (x) 


1 

v/?ir 

1 

VT* 


/; 

£ 


♦ + (e)e 


-16X 


♦ (B)e 


-iBx 


ds. 


de. 


(43) 


(44) 
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Note that the Integral path Is here the real axis, which is In the region spe- 
cified by equation (42) 

The Fourier transform of equation (35) yields 

♦ + (b) ♦ b_(b) - -Ika a r m ( B ^[ F mn (a) + M®)]* ( 45 ) 

Here r m , which is obtained in appendix A, Is the Fourier transform of Ym : 

r " ( ' ,) (46) 

And 






l 

V5V 



dx ^.(a, x)e^ Bx . 


(47) 


The integral in equation (47) is defined only for 


*•(') > - • 


However 


F mn ^fincd ’ n ent i re B space except at the pole b 
stitutes the analytic continuation of the integral for 9m( b) 
Rearranging equation (45), one obtains 



, and con- 



L(8)[» + (b) ♦ F mn (e)] - -Mb) ♦ F mn (e). (48) 


where 


L(b) * 1 ♦ 1VZ7 ka a r m ( a , a ; b ) . (49) 

As shown in appendix A, r m (a,a;a) — 1 — . 0, thus 

L(b) |b I --I 1. (50) 

The next step of the Wiener-Hopf technique is to decompose equation (48) 
into two factors, one analytic in the UHP and the other analytic in the LHP. 

To this end, L(b) shall be factorized first. With the substitution of equation 
(All), equation (49) is written as 

aJ'(a) - 1kan 2 J_(a) 

L(e > - .j « ■ «,&.) • (51 » 


where 
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The numerator and the denominator In equation (51) are even functions of a 
for m even, and are odd functions of a for m odd. Thus, for m even the 
two factors are entire functions of 0 . For m odd, the two factors are con- 
currently multiplied by a to be made entire functions of 0 . In either case, 
th^ ^ntire functlns possess simple zeroes and equation (51) can then be written 


L(B) 



;nr 

mi J 


1 ♦ 


mi J 


where 


, J ; (ka > - Vm< ka > 

L 0 ' Jftlca) - UjJJka)' 


The factorization of L(b) immediately follows: 

L(e) - L + (b)/l - (0) 

where 


L + (0) 



(52) 


(53) 


(54) 


(55) 


L_(0) 


D f ■ # ■ 


mt J 


(56) 


Note that 

L_(-B) - Lq/L + (b). 
Equation (48) can now be written as 

M*> [»♦<»> * F mn(«)] - L - Wn U ) F m <»> 


(57) 


L -(»)[-‘-C) * F mn(»)] - L -( k ») W" < M ) 
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Here the term L_ (-k^) • F^e) been subtracted from both sides to re- 
move the pole of F^fa) from the right side. 

The left and right sides of equation (58) are analytic, respectively, In 
the UHP and the LHP. The two sides are the same and analytic In the common 
analytic region specified In equation (42). Thus they are analytic continua- 
tions to each other,! 4 ) and equation (58) defines a runcllon which Is analytic 
in the entire a space. As a ♦ •, L*(a) ♦ constant and «*(a) ♦ (1/a). 

The asymptotic behavior of k(b) Is Inferred from the behavior of v(x) near 
x - 0.'°) It follows then that both sides of equation (58) tend to be zero 
like (1/e), as a ♦ -. According to the Liouvllle's theorem,' 9 ) the function 
must equal zero everywhere, and one obtains 


♦.(«»■ [ L - 

(59) 

•-<»- -[ L - (- k ii ) )/ L » - 1 ] F „ ( *)- 

(60) 


Reflected and Transmitted Waves 

The reflected wave is determined from equation (44) with the substitution 
of equation (60). The integral in equation (44) can be replaced by a contour 
integral along a path C closed in the UHP as illustrated in figure 3: 

» (x) « — <f ♦ (e)e' ie * da for x < 0. (61) 

y/li * / C 

The function *_(e) is singular only at the simple poles a - k^ j - 0,1,2, 
which corresponded to simple zeroes of L_(a). Thus, equation (61) involves a 
Cauchy type integral, which is readily evaluated as 


r_(x) 



for x < 0, 


(62) 


where 



(a; 
y mn v ' 

y mj vay 




( 63 ) 
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Equation (63) determines the reflection coefficients, which have been defined 
in equation (17). 

The transmitted wave is determined from equation (43) with the substitu- 
tion r c equation (59). The integral in equation (43) is replaced by a contour 
integral along a path C' shown in figure 3. The function * + (s) is singular 

( 2 ) 

only at the simple poles at e « “ k mv * v ■ 0,1,2 corresponding to the simple 
zeroes of L+(p). One can then readily obtain 


v + (x) 



i(2) x 

O 1 * 2 >C’ < a > e m 




for x > 0, (64) 


where 


T n U 

nv ' 


- 2 ) 


oli’u) 

32 ^ 


— n — 


k ( 2) 

mv 

k (») - k £) • 
mn mv 


(65) 



The second term in the right side of equation (64) is none other than the incident 
wave, thus for x > 0 the net field is, from equation (29), 


v(a,x) = 



ik < 2 >* 


T n <* * a)Q* 2 ) (a)e "* 

nv mv 


for x > 0. 


( 66 ) 


This is the transmitted wave, and equation (65) determines the transmission 
coefficient, which has been defined in equation (18). 


DISCUSSION: RECIPROCITY AND LIMITING CASES 

Foregoing detailed numerical calculations, we will discuss here the lim- 
iting cases of a small admittance shift a and of I9W frequencies. The dis- 
cussion will also include the reciprocity relations.^) 


Reciprocity Properties 

As in appendix B, the infinite products that are explicitly comprised in 
equations (63) and (C' can be expressed in terms of the factor L+ which can 
be evaluated numerically.^) The reflection and the transmission coefficients 
can then be written as 
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ak, 


HT 

mj 


-ik a 
TJTj 
mn 


+ k, 


rn 

mj 


L udLd) 
L O u mj u mn 


rwsT 


(67) 


O 1 * 21 - 


ak, 


m 

mv 


ik a 

?*T 

mv 


L + k 


- k 


nr 

mn 


mv 


( 1 ) ( 2 ) 

U U 

*mn M mv 

tit; 

mn / 


(68) 


where we have also used equation (57). 

The reciprocity relations can be readily established from equations (67) 
and (68): 


kiVR”.(D - k£i , n7_(i>. 


.<U R " 

'mj nj' 


mn jn' 


(69) 


l<(2)ym 

mv nv 


(1 


2) = k 


(l) T m 
mn vn 


(2 


1 ). 


(70) 


Here T™ n (2 -► 1) is the transmission coefficient for the (m,n) mode generated 

in the duct 1 upon the incidence of the (m,v) mode from duct 2. Recall that 
reciprocity is a characteristic of a system linear vibration. For a more de- 
tailed discussion, one may refer to reference 10. 


Small Admittance Shift 

When the admittance shift is small such that I ka a l « 1, equations (67) 
and (68) become simple expressions: As detailed in appendix C, the reflection 

and the transmission coefficients are written as 



(71) 



2 ) * 1 - 


ik a 

TaT 


"u (1) 

M mn 

rnr 

_ mn J 


2 


+ 0 


[(ka a) 2 ] 


for v = n. 


(72a) 


ik 


* 


ak 


(1) 17TT7 : k '(17 


+ o 


mv 


mv 


[(ka a) 2 ] 


for 


4 n. 


(72b) 


mn 


As expected, a very small admittance shift causes very little reflection or 
very little mode scattering in the transmission. As a increases with a fixed 
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ni, the reflection coefficient increases linearly with a. The behavior of the 
transmission coefficient may be discussed in two separate cases. For v = n, 

the transmission coefficient deviates linearly with a from unity. For v * n, 
the transmission coefficient linearly increases with a, and is inversely pro- 
portional to - kji^. 


Low Frequency Limit: ka << 1 or I a I 


In this case, one may be interested in the incidence of the fundamental 
mode only. For simplicity, we also assume = 0, that is, the duct element 1 

is hard-walled. The reflection and the transmission coefficients are then ap- 
proximated as 

“oo (1) ' ' [ x ' fr! * 0 (r)]> (73) 

T So (1 * 2 > s yrr*°(!r)- (74) 

As ka + 0, the wave will be completely reflected from the impedance discon- 
tinuity. As the frequency increases from zero, the amplitude of the reflection 

coefficient decreases from unity as VlTa. It is interesting to compare this 
with the low frequency behavior of the reflection from the open end of the 
duct. For the latter case,(^) the reflection coefficient is 




(75) 


CONCLUDING REMARKS 

The Wiener-Hopf technique has been employed to study modal scattering of 
sound in a circular cylinder subjected to an axial discontinuity of the wall 
impedance. Analytical expressions have been derived for the reflection and 
the transmission coefficients. The analytical expressions display some symme- 
try properties such as the reciprocity relations, and can be immediately eval- 
uated by numerical integration of the factors in the reflection and the trans- 
mission coefficients without calculating the complex modal eigenvalues. 

The quantitative discussion was confined to the limiting cases of small 
admittance shift, and of low frequency. With the small admittance shift, the 
reflection and the transmission coefficients depend linearly on the admittance 
shift . 

The case of low frequency includes the fundamental mode that is incident 
from a hard-wall duct and reflected at the junction with a soft-wall duct. At 
the very low frequency, the reflection coefficient is almost unity in the am- 
plitude and with phase charge of 180 degrees. This result is similar to the 
low frequency limit for the reflection of the fundamental mode from an open 
end of a hard-wall duct. However, the deviation from unity with increasing 
frequency is quite different for the two cases. For the present problem, the 
deviation is proportional to the square root of the frequency, whereas, for 
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the reflection from the open end, the deviation is proportional to the square 
of the frequency. 

Another low frequency characteristic is that the deviation of the reflec- 
tion coefficient from unity is inversely proportional to the square root of 
the admittance of the soft wall. 
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APPENDIX A 

FOURIER TRANSFORM OF GREEN'S FUNCTION 

The function r m can be obtained directly from equation (28) by using 
equation (46). The procedure Is, however, cumbersome. We employ a more ele- 
gant method In which r m Is obtained as a solution of the Fourier transformed 
equation. The equation for the reduced Green's function Is 


l r ar ar 


m 2 » 2 .2 

T ~ 7 k 

r ax 


A , «(r - r„)« 

I r m (r.r 0 ;x - «„) - 


(x - X n ) 


Consider the integration 


(Al) 


- - — f eq. (Al)e^ C d; 

\/Ti -« 

where c m x - xq. Effecting this integration, one obtains 
^1 d 




id M d nt i o 1 » \ 

F3r r 3F'T t) r m (r * r 0 ;e) " 
r 


2 2 \ s(r - r Q ) 


\/2* 


where 


a 2 - (k 2 - e 2 )a 2 . 


The solution can be expressed as 


r m (r,r 0 ;e) 


for r o < r 


V?* r o W ^ y l’ y 2^ 

r*i 


o y 1 (r)y 2 (r 0 ) 


for r < r. 


(A2) 


(A3) 


(A4) 


(A5) 


Here yj and y 2 are homogeneous solutions of equation (A3), and W is the 

Wronskian. It follows from equation (A5) that y,(r) should be analytic at 
r - 0, and thus one has 1 


n< r > ■ (» j) • <a«) 

The solution \/2 linearly independent of yj, and is here chosen as 

>2 ' J m (« 0 * 8 »„ (« l) ' 

The constant B Is determined from the boundary condition at r ■ a as given 
in equation (13), which was imposed on the Green's function: 
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dy 2 

- jy.— - ikn^ * 0 at r * a. 

On inserting equation (A7) into equation (A8), one obtains 

B • - «y;u) - ik.njVj.j- 
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m 


The Wronskian is obtained as 


W 


2B_ 

" r 0 ' 


Summarizing equations (A 5) to (A10), one obtains for r = rg = a 


m 


r m ( a » a * e ) o J^(a) - IkanjJJa)* 


(A9) 


(A10) 


(All) 


Using the asymptotic form of Jm for a large argument; (12) one can show 

r n| (a.a;s)l£L-lo(i). (A12) 


18 



ORIGINAL PA.lv.' IH 
of poor quality 


APPENDIX B 

INFINITE PRODUCT TERMS 


The infinite products contained explicitly in equations (63) and (65) will 

be expressed in terms of the factor L+. 

Let Dj denote the infinite product in equation (63): 




r kun 

1 ■ ziff 

] 



" k dn 




1 

t4i 


mt _ 


(Bl) 


Compared with equation (56), this product can be written as 


1 im 


1 


*1 " k (1 > 

* * K mj " 




6 

T 

mj J 


(B2) 


Using equation (54) and L+(e) being analytic in the UHP, one can write equa- 
tion (B2) as 


D, 


1 


1 im 


L + k , 


0 - k 


(1) 

mj 


L(B) 


1 - 


B 

T 

mj J 


k^V 


(B3) 


As noticed from equation (51) or (52), L(b) possesses a simple pole at 
B * k mj^ ‘ Near the po ^ e » *-(b) can be ex P an ded as 


1(B) 


k (i) 

mj 

' (If 

. M mj . 

z 

J' 

mj m 

MS’) 

- ikan 2 J r 

. (45 1 )] 



mj m \ 

w 

B - 

mj 



(B4) 


On inserting equation (B4) into equation (B3), one obtains 

-|2 


ika a 


J 1 




(i) 


mj 


L aK mj J 


(B5) 


where we have used 

4^ (4i>) - (4j>) - 

Consider now the infinite product in equation (65): 


(B6) 
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SMALL ADMITTANCE EXPRESSIONS 

Here we will determine the limiting values of the eigenvalues, the 
propagation constants, and the factor L+. These limits have been used to 
derive equations (71) to (74) from equations (67 and (63). 


Small Impedance Jump: |ka a I « 1. 

Consider the eigenvalue equations derivable from equations (13) and (14): 


“ (1) J m (« (U ) - ika "l J », 1 

;° (1) ) . o. 

(Cl) 

(' ,l2) ) - ,k Vm 1 

(« (2) ) - o- 

(C2) 

Set 

(2) (1) c 
a - a = e • 


(C3) 

In the limit of |ka a 1 << 1, one expects that 

|e|« .<»|. 


(C4) 


To determine e, equation (C2) is expanded around and equation (Cl) is 
used to obtain an algebraic equation for t. With the first order expansion, 
one obtains 


-ika a 



c mi * nr 

a mt 


(C5) 


where the mode numbers (m,t) are used as the subscripts on e and a(*). 
Equation ( C 5) obviously satisfies the inequality (C4). The equation is not 

valid for the case of = 0 which is separately treated later. Also it 

should be mentioned that the derivation of equation (72a) requires the second 
order expansion, which will not be discussed here. 

For the propagation constants, one sets 


a 




mf 


(C6) 


Using equations (10), (C3), and ( C 5) , one obtains 
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■'“■“Sro 


(C7) 


mi 


Unless the frequency Is close to the cutoff frequency of the (m,i) mode, 
equation (C7) satisfies the inequality 


M 


« a k. 


,(D 

'mi 


Set 


where 


" L 0 9 2 (B)/g 1 (6) 


g s 00 


OE 


+ » 
mi J 


s « 1,2. 


Using equation (C6), one may obtain 


9 2 (b) * 


t*0 




1 ♦ B 
1 7TT> 


bk 


mi 


mi a k 


* 9 1 (B) <1 - 



3T F 5^} 


With the substitution of equation (C7), equation (Cll) is written as 

9 2 (b) » g x (B) [l - ika a x ( b )] , 

where 


x(b) 



uvr 

7 

k df 

w mi 


1 + JB1_ 

B 

TOT 

/ 

_mi_ 

/ 

_ 


(C8) 

(C9) 


(CIO) 


(Cll) 


(C12) 


(C13) 


The function x(b) is a finite and well behaved function of e in the UHP. 
For equation (53), 
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L 0 ■ 1 - araykaJ^MfanjJ^k*) - (C14) 

From equations (C9), (C12), and (C14), 

L + (a) - 1 - ika a [x(e) - ra;TkT) "- 1 < E a^'J "( Va )J * (C15) 

On Inserting equation (C15) Into equations (67) and (68), one obtains equations 
(71) and (72b). As mentioned earlier, the second order expansion Is required 
for equation (72a). For equation (72b), we have also used 

“ra! 1 * “ii 1 [ l * °< kJ *)] • < C16 > 


Low Frequency Limit: ka « 1 and lal 


As in the text, we consider the case of nj = 0 and of the Incidence 
of the (o,o) mode (the plane wave mode). We have then 


a » n 2» <*00^ = °» k oo^ 


k. 


Set <*^ * c, then the expansion of equation (C2) yields 


and 


g * V-2ika & 


ak^ * y/ZWTT, 


(C17) 


(C18) 

(C19) 


i - 2in 

L 0 a YT> 


(C20) 



L 


+ 



* 2 



( 2 ) 

p 00 




1. 


(C21) 


(C22) 


(C23) 
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DUCT ELEMENT 1 — J— DUCT ELEMENT 2 



pcu^p*^' * L pcu r /p-J?2 


Figure 1. - Circular cylinder with admittance shift 
•t x - Ct im.nlmcde incident; im.ji modes reflected; 
)m,uf modes transmitted. 




